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Abstract. In a previous work we apply lattice point theorems on hyperbolic spaces obtaining 
asymptotic formulas for the number of integral representations of negative integers by quadratic 
and hermitian forms of signature (n, 1) lying in Euclidean balls of increasing radius. This 
formula involves an error term that depends on the first nonzero eigenvalue of the Laplace- 
Beltrami operator on the corresponding congruence hyperbolic manifolds. The aim of this 
paper is to compare the error term obtained by experimental computations, with the error 
term mentioned above, for several choices of quadratic and hermitian forms. 

Our numerical results give evidences of the existence of exceptional eigenvalues for some 
arithmetic subgroups of SU(3, 1), SU(4, 1) and SU(5, 1), thus they contradicts the generalized 
Selberg (and Ramanujan) conjecture in these cases. Furthermore, for several arithmetic sub- 
groups of S0(4, 1), S0(6, 1), S0(8, 1) and SU(2, 1), there are evidences of a lower bound on the 
first nonzero eigenvalue better than the already known lower bound for congruences subgroups. 



1. Introduction 

Let Hp denote the F-hyperbolic space (for F = M, C or H) of dimension n. Let x, ?/ G H|J 
and let F be a discrete subgroup of the isometry group of HJJ such that r\HJJ has finite volume. 
Hyperbolic lattice point theorems give asymptotic formulas (as r — )■ +00) for 

N.,y{r) := #{7 e r : rf(x, 7 ■ < r}. 

This problem has been investigated by several authors (see for instance |Pa75j . |LP82] . |Le87j . 
|PR94j . |BMW99] . |GN12j ). most of them concentrate on real hyperbolic spaces. The main 
theme in these works links the asymptotic behavior of Nx^y{r) with spectral data of the Laplace- 
Beltrami operator on r\H5J. 

For example, if F\H|[J is not compact, Bruggeman, Miatello and Wallach |BMW99"] showed 
that 

(1-1) N^Ar) = XJoKT\^) + ^^-^"^^^-^^^ '"'"''^ + ^ (e(^^^+^)^) 

for all e > 0, where 

r = rF = dimM(F) = 1,2,4, 
P = P„,F = ^r-l = i^,n,2n + l 

for F = M, C, EI respectively, vol(S'™) denotes the volume of the m-dimensional unit sphere, 
Bi, . . . , i?Ar are positive (explicit) coefficients, the numbers 

< Ai < ■•• < Atv < 
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are the exceptional eigenvalues of the Laplace-Beltrami operator on r\Hp counted with multi- 
phcities and, for each 1 < j < A^, 

and ipj is the normahzed eigenfunction associated to Xj. We recall that f{x) = 0{g{x)) if and 
only if there are M,C > such that < C g{x) for all x > M; this is also denoted by 



f{x)^g{x). If F = M, (1.1) holds already with e = by |Le87] . In general, A-,- and (pj are not 



determined for any j, hence it is customary to work with the abbreviated formula 

2i-"vol(S™^ 
2pvol(r\H|^ 



(1-2) N.Ar) = KZlTZ'n? e^^- + Q(e-) 



where 

(1.3) r = T{T):- 




ifz.i+p>2p^, 
otherwise. 



We can see that the exceptional eigenvalues play an important role in the distribution of 
the lattice points. The existence or not of exceptional eigenvalues is an open problem. By 
definition, Ai(r\Hp) > p^p if and only if r\Hp has not exceptional spectrum. The Selberg 
eigenvalue conjecture says that 

Ai(r\H2)>p|j, = l/4 

for all congruence subgroups F in SL(2, Z). This is part of the Ramanujan-Petersson conjecture 
at infinity. Selberg |Se65] established the bound Ai > 3/16 by using the Kloosterman-Selberg 
zeta function and Weil estimates for Kloosterman sums. The best lower bound available at the 
moment is Ai > 975/4096 ~ 0.238 . . . proved by Kim and Sarnak [KS03J . Furthermore, the con- 
jecture has been verified in some explicit cases: Roelcke |Ro56j for F = SL(2, Z) (see also |DI83j ). 
Huxley |Hux86j for Fo(A^) for all < 19, and more recently, Booker and Strombergsson |BS07j 
for Fo(A^) for all < 857 squarefree. 

It is also expected that this conjecture holds for n = 3; in this case the conjecture says 
that Ai(F\H|) > pIk = 1 and Sarnak |Sa83] proved that Ai(F\H|) > 3/4 for all congruence 
subgroups of SL(2, (9), where O is the ring of integers of an imaginary quadratic number field. 
Some particular cases are verified in |EGM89j . 

It is known that the generalization of Selberg's conjecture in higher dimensions is false. When 
F = M and n is arbitrary, for Q a quadratic form with rational coefficients of signature (n, 1), 
we say that a subgroup of 

(1.4) Fj := GL(n+ 1,Z) nSO°(Q) 

(see Subsection 2.1 for arbitrary F) is a congruence subgroup if it contains Fq(A^) := {g G 
T^Q : g = In+i (mod A^)} for some A^ G N (see Definition 2.1 for arbitrary F). One can find in 
|CLPS9lt §6] an explicit construction via theta-lifting from SL(2) of exceptional spectrum of 
Fq(A^)\H^ for large A^. However, Selberg and Sarnak's lower bound were generalized to 

9r7 — 3 

(1-5) Ai(f\h;^) > 

for all n > 3 and for any congruence subgroup F. This was proved independently by Elstrodt, 
Grunewald and Mennicke |EGM90j and Cogdell, Li, Piatetski-Shapiro and Sarnak |GLPS9l] . 



In the unitary case the situation becomes even darker. Jian-Shu Li |Li91j (see Theorem 2.6) 
proved that the first eigenvalue Ai for the Laplace-Beltrami operator on the subspace of non- 
degenerate forms (nontrivial Fourier coefficients) in L^(F\H£) satisfies Ai > 2n — 1, where F is 
any congruence subgroup of SU(n, 1). 
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The aim of this paper is to check ( 1.1 ) numerically by assuming that there are no exceptional 
eigenvalues in F^yHp for several choices of Q. Our computations seem to contradict this 
assumption (Ai > p^) for some hermitian forms Q of signature (n, 1) with = 3,4, 5. In other 



words, we have numerical evidences to affirm that F^yHp has exceptional spectrum. 

We also consider quadratic forms of signature (n, 1) for n = 2,4, 6, 8 and hermitian forms of 
signature (2, 1). In these cases, our numerical results do not seem to contradict the assumption 
Ai > Hence, our calculations give evidences of the weaker condition 



:i.6) Ai(rJ\H^) > 



An 



(n + 1) 



2 

2 Pn,F 



if and only if p + z/i < 2p 



n 



n 



1 



since, by (1.3), this is the minimum lower bound that ensure the smallest error term in (1.2). 



If (1.6) holds (e.g. there are no exceptional eigenvalues), then we have that (1.2) holds with 



:i.7) 



r 



2p 



n 



n + 1 



+ e 



n{n — l)/(n + 1) 

2n^/{n + l)+e 

{An + 2)n/{n + 1) + e 



if F 
if F 
if F 



C, 



for every e > (recall that £ = when F = M). To check ( |1.1 ) in a numerical way, we use a 
recent arithmetic application. 

Hyperbolic lattice point theorems have been applied in different contexts, mostly on arith- 
metic problems (see |EGM88j ). For instance, Ratcliffe and Tschantz |RT97j . by applying the 
results of [ LP8 2] and |Le87] . obtained the asymptotic formula, as r — )■ +oo, for the number of 
representations of a negative integer by the Lorentzian quadratic form that are contained in the 
ball of radius r centered at the origin in M*^"*"^. In |Lal2j . we generalized jRT97j to more general 
integral quadratic forms of signature (n, 1) and also to hermitian forms. We next introduce 
some notation to state the main formula of |Lal2j . 

For F = M or C, let O denote a maximal order in F, that is, (9 = Z in the real case and 
the ring of integers of an imaginary quadratic extension of the rational numbers if F = C. We 
consider, for n > 2, an integral F-hermitian matrix {B G M(m, O) such that B = B*) 



Q 



A 



—a 



with a G N and A G M(n, O) a positive definite integral F-hermitian matrix. We also denote 
by Q the induced F-hermitian form Q[x] := x*Qx = A[x] — a\xn+i\'^ {x G O"''^^ and x : = 
{xi, . . . ,x„)*). For t > 0, put 

(1.9) Nt{Q, -k) = {xE : Q[x] = k, \xn+i\ < t}. 



In jLal2j . by applying the formula (1.1) from |BMW99] to the discrete subgroup Tq given by 
(1.4), we establish the following formula. 



Theorem 1.1. If —k is represented by Q, then 
(1.10) Nt{Q,-k) = C{Q,-k)t^^ + 0{t^) 

where 

2(r-l)(n+l) V0l(5"''-1) W2 



as t ^ +00, 



:i.ll) 



C{Q,-k) 



MoT^^ I det Q|i 



F( 



6{Q,-k). 



Here, do denotes the discriminant of the quotient field ofO, 6{Q,—k) = Ylp^piQ^ 



local density of the representation Q[x] = —k (see (3.1)) and t is as in (|1.3). 
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Now we can explain our strategy in more detail. We set 

Nt{Q,-k) 



[1.12) 



C{Q,-k) 



thus ^(t) ^ t'^^'^P and limt_>oo ^(^) = since r < 2p. With the help of a computer, we 
quantify \E'(t) for large t, estimating its decay as \l/(t) ~ Bt"^ (i.e. limj_>oo '^{t)/t'^ = B) for some 
cr < and 5 > 0. Next, we compare the obtained error term 0{V~'^'') of '^{t) in (1.10) with 
our estimate 0{t''). As we mentioned before, we are mostly interested in comparing them by 
assuming that Ai(rQ\Hp) > pn,¥ (i-e. there are no exceptional eigenvalues in F ^yH g), or more 
generally by assuming the weaker condition (1.6); in both cases, r is given by (1.7). 

We state the conjectures evidenced from the numerical results. Let 
identity matrix. 



In denote the n x n- 



Conjecture 1.2. Let Q be the C-hermitian form diag(J„, — a) over for n = 3,4,5, 

where a G Z and 1 < a < 15. Let Tq C SU((5, C) = SU(n, 1) be the associated discrete subgroup 
given by (1.4) and let Ai be the first nonzero eigenvalue of the Laplace- Beltrami operator on 
TqXH^. Then Tq\H.^ has exceptional spectrum (i.e. Xi < p'^^ = n^). Moreover, 

Conjecture 1.3. Let Q be the ¥-hermitian form diag(/„, — a) for ri = 4, 6, 8 z/ F = M and, 
over Z[a/^] for n = 2 if¥ = C, where a e Z and 1 < a < 15. Let Tq and Ai as above. Then 

A > 

where pn,¥ = ^ for F = M, C respectively. 

We have made similar computations for other quadratic forms of signature (4, 1) and (6, 1). 



2 

\2 Pn,¥-i 



They are listed at the end of the paper. Conjecture |1.3| also holds for these cases. 

This paper is organized as follows. In Section |2| since Fq is trivially a congruence subgroup 
(see Definition 2.1), we obtain upper bounds for r by applying known lower bounds for Ai 



for general congruence subgroups. Sections [3] is devoted to compute, for arbitrary diagonal 
nondegenerate F-hermitian forms Q and any G Z, the local density 6{Q,k) = Ylp^piQ^ 
In Section |4] we compute the main coefficient C{Q, —k) and also we give an approximation for 
the number a. In Section [5] we present tables with the numerical results and their respective 
conclusions. All the necessary computations have been done in |Sage| . 
Let us conclude this introduction by looking at two illustrative examples. 



'4,1 



( ^'^ _j \ where 



Example 1.4. Consider the M-hermitian form, or 'quadratic form', Q 

I4 denotes the 4x4 identity matrix. Here we have C{Q, —1) = 5 (see Subsection 4.1). Fur 
thermore, (|L5|) says that Ai > 5/4, thus r = i^i + p4,R < ^(3/2)2 - 5/4 + 3/2 = 5/2 from ( [L3| ) 
(see also Corollary 2.4). Hence, we have 

Nt{h,i,-1) = 5t^ + 0{t^-^) ast^+00. 

On the other hand, if there do not exist exceptional eigenvalues of the Laplace-Beltrami operator 



on Fq\H^ (or more generally (1.6) holds), (1.7) ensures that 

Nt{h,i,-1) = 5t^ + 0{t'^-^) ast-^+00. 

The formulas given above mean that \l/(t) ^ and \l/(t) ^ respectively. On the 
other hand, our approximation of a is —0.93, more precisely, ^'(t) ~ 13t~°-^^. Fi gure 1 shows 
the points (m, \l/(m)) up to 10^ and the three different curves mentioned above. 
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Figure 1. Graphic of "^tiQ, -1) when F = M, Q = h,i and k = 



1. 



5 






The graphics include the points (m, \l/(m)) for m G [10^, 10^] above and for m G |10^, 10^] 
below, curves of decreasing order t~^/^ and The curves are chosen so that they 

have the same value at the endpoint on the left. 



Pi 



We can see that the approximating curve is under the curve obtained by assuming Ai > 
]g = 9/4, thus these computations do not contradict the assumption of the nonexistence of 
exceptional eigenvalues in rQ\Hj|. On the other hand, this numerical test gives an evidence for 
the lower bound Ai(ro \H^) > pl^An/{n + l)^ = 36/25 = 1.44. 



Example 1.5. We now consider, for F = C and O 

18 (see Subsection 



We have C{H,-1) 

Nt{H, 



Z[-\/— 3], the hermitian form H = I^^i 



4.1 ), hence 
1) = 18ti" + 0(r) ast^+cx). 
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Figure 2. Graphic of ^t(g, -1) when F = C, C = H = h^i and A; = 1 





The graphics include the points (m, \I/(m)) for m G |10^, 10'*] above and for m G |10^, 10^] 
below, curves of decreasing order and The curves are chosen so that they 

have the same value at the endpoint on the left. 



where r is given by (1.3) 



In this case, if we assume that there are no exceptional eigenvalues, or more generally (1.6), 
then the above formula holds with r = 50/6 + e (50/6 ~ 8.33) for each e > 0, hence \E'(t) -C 
^-i.67_ However, this assumption contradicts our approximation \E'(t) ~ 56 thus this 
is numerical evidence on the existence of exceptional spectrum in r5^\H|.. Moreover, the 



estimation r — 2p ^ a = —1.36 tell us that Ai ~ 11.82 since r = v5^ — Ai + 5 and 2p = 10. 

Figure [2] shows the points (m, \l/(m)) up to 10^ and the three different curves mentioned 
above. 



A NUMERICAL STUDY ON EXCEPTIONAL EIGENVALUES 



7 



In this case, Jian-Shu Li's lower bound |Li91j (see Theorem 2.6) imphes that Ai > 9, where 
Ai denotes the first nonzero eigenvalue for the Laplace-Beltrami operator on the subspace of 
nondegenerate forms (nontrivial Fourier coefficients) in L^(r^\HJ^). If we assume that this lower 



bound holds on all L2(r^\H^), we obtain that r < 9 (see Remark [2^, hence \E'(t) G 0{t~'^] 
Therefore, our estimation \E'(t) ~ 56 give evidence on the extended Li's lower bound for 



^ H- 



2. Upper bounds for r 



2.1. Hyperbolic spaces. In this subsection we may assume that F = M, C or H. It contains a 
brief review of Sections 2 and 3 in |Lal2] . We consider F""*"^ as a right module over F. Let Q be 



an F-hermitian matrix as in (L8). This matrix induces the F-sesquilinear form Q(x, y) := x*Qy 
and the F-hermitian form Q[x\ := Q{x,x) = x*Qx. 

We consider the n-dimensional F-hyperbolic space realized by the Q-Kleinian model 

(2.1) H^(Q) = {[x] gPF":Q[x] <0}. 

Under a standard Riemannian structure described in |Lal2[ §2], the distance between two points 
[x] and [y] in Hp((5) is given by 

cosh(rf([x],M)) = ^J^j|^^^ 

Let U(Q,F) be the Q-unitary group and let SU((5,F) be the special Q-unitary group, that is 

U(g, ¥)={ge GL(n + 1, F) : Q[g] = g*Qg = Q}, 

SU(g,F)={(7GU(g,F):det((7) = l}. 

The action of these groups on Hp(g) is given by g ■ [x] = [gx] for g G U(g,F). This action is 
by isometrics, transitive and the only elements that induces the identity map are those in the 
center Z(U(Q,F)) of U(Q,F). 

Let G = SU°(g,F) be the identity component of SU(Q,F). When F = C or H, the group 
SU(Q,F) is already connected, thus G = S\J{Q,¥). However, if F = M we have 

G = PSO(g) := {g G U(g,M) : (7n+i,n+i > 0}. 

In all cases, Hp(Q) coincides with the symmetric space G/K, where K is the isotropy subgroup 
in G of [e„+i] G H^(Q) (see [LiT2l §3]). 

Let Tq = U(Q,F) n M(?7, + 1, O) denote the unimodular matrices in U((5,F), and let Fg = 
Fq n G. It is a simple matter to see that 




c, 



where u is the number of units in O. These groups act discontinuously on Hp((5) and the 
quotients Fq\H|^(Q) and Tq\H.^{Q) have finite volume. Moreover, Tq\B.^{Q) is compact if 
and only if the F-hermitian form Q is anisotropic over the quotient field of O. 

2.2. Exceptional eigenvalues. Let A be the Laplace-Beltrami operator on }i^{Q). For a 
general discrete non-cocompact subgroup of finite covolume of G, the operator 

-A : L\T\}l-m L\T\ll-m 

is known to be essentially self-adjoint and positive and one has the spectral decomposition 

L\T\E^m = Ll{T\}l^m L^(F\H;^(Q)), 
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where the spectrum of A is discrete (resp. continuous) in L^(T\H'^{Q)) (resp. L^(r\Hp(Q))). 
Let p = (n + l)r/2 — 1 where r = dim^{¥). There are only finitely many eigenvalues of the 
discrete spectrum of A (counted with multiplicities) such that 

< Ai < ■ ■ ■ < < pi 

These are called exceptional eigenvalues. 

Definition 2.1. Let Q be an F-hermitian form with coefficients in the quotient field of O with 
signature {n, 1). A subgroup of Tq is called a congruence subgroup if it contains 

r° (a) := G : ^7 = /n+i (mod a)} 

for some ideal a of O. 



We recall that the error term in the main formula of Theorem 1 . 1 
is given by (L3) associated to the discrete subgroup Tq. Hence, a 



is 0{V) where r = r(ro ) 
ower bound for Ai makes 



the error term smaller. Our goal in this section is to apply known lower bound for Ai(r\Hj[J) 
with r a congruence subgroup, obtaining upper bound for r. 

We begin in the hyperbolic plane case (i.e. F = M and n = 2). The following result is due to 
Kim and Sarnak |KS03] . 

Theorem 2.2. For any congruence subgroup T < SL(2,Z), the first nonzero eigenvalue Ai for 
the Laplace- Beltrami operator on r\H^ satisfies Ai > 975/4096 ~ 0.238 .... 

In Introduction we mentioned that, for any n > 4, there are congruence subgroups such 
that their associated hyperbolic manifolds have exceptional spectrum (see |CLPS9l ', §6]). The 
following result is the optimal lower bound for n > 4. It was proved independently by El- 
strodt, Grunewald and Mennicke jEGM90[ Thm. A] and Cogdell, Li, Piatetski-Shapiro and 
Sarnak [CLPSM]. 

Theorem 2.3. Let n > 3 and let Q be a quadratic form with rational coefficients such that 
Q is of signature (n, 1) and isotropic over Q. For any congruence subgroup T < Tq, the first 
nonzero eigenvalue Ai of the Laplace-Beltrami operator on r\H]^((5) satisfies 

2n-3 
A.> — 

Corollary 2.4. Let F = M, let T be a congruence subgroup acting on the n-dimensional real 
hyperbolic space, and let t = T(r) given by (1.3) associated to T. Then 

n — 3/2 if n > 3, 



(2.2) 

Proof. We recall that Uj :- 



T < 



2/3 



ifn = 2. 



— \j. For n > 3, Theorem 



2.3 



1^1 < 




2n 



n 



imply that 
2 



4 



From (1.3) we conclude that t < ui + p < {n — l)/2 + {n 
For n = 2, Theorem 2.2 implies that ( 1.6[ ) holds since 

4np2 _ 2 

(n + 1)2 ~ 9 

then r = 2pn/{n + 1) = 2/3 by (O). 



2 

2)/2 



n 



3/2. 



975 

0.22- ■■ < < Ai, 

4096 - ' 



□ 
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Remark 2.5. Note that Corollary 2.4 gives the best error term for (1.10) that can be obtained 
from jLal2] . i.e. having used the lattice point theorem of Levitan jLe87j . 

Bounds of this kind are not known in the complex case. A related result was proved by 
Li jEoTl Cor. 1.4]. 

Theorem 2.6. Let n >2 and let T be a non-cocompact congruence subgroup ofSl]{n, 1). The 
first nonzero eigenvalue Ai for the Laplace- Beltrami operator on the space of nondegenerate 
forms in L^(r\H^(Q)) satisfies 

Ai >2n-l. 



Remark 2.7. The analogous procedure to that in Corollary 2.4 in this case, by assuming that 
Theorem 2.6 holds for every form in L^(r\H£((5)), gives that vi < n — 1, hence 



(2.3) T<vi + p = n-l + n = 2n-l. 

3. Calculation of local densities 

The purpose in this section is to compute the local density 6{Q, —k) of the representation 
Q[x] = —k. Here we assume F = M or C. Recall that r = dimK(F), (9 is Z in the real case 
and the ring of integers of an imaginary quadratic number field if F = C, and do denotes the 
discriminant of the quotient field of O. 

For an arbitrary nondegenerate F-hermitian forms and any /c G Z, the local density of the 
representation Q[x] = k is defined by 

(3.1) S{Q,k) = l[ 6p{Q,k), 

p 

where the product runs over all positive prime numbers and 

, , # G (O/p^Oy" : Q[x] = k modp^] 

Throughout this section, we only consider diagonal integral nondegenerate F-hermitian forms 
of rank m, thus, they are given by matrices of the form 

(3.3) I •.. J (a^ gZ, a^^OVz). 

\ am/ 

Note that the signature of this matrix is (mi,m2) where mi (resp. denotes the number of 
positive (resp. negative) entries in the diagonal. 

As usual, let ({s) denote the Riemann zeta function. We shall use the well known formula 

3.4 CM = ^"n — ' 

2 ml 

for every m odd integer, where B^. denotes the m-th Bernoulli number 

We will denote by (— ) the Kronecker symbol. If Z) is a quadratic discriminant (i.e. D = 1 
(mod 4) square free or D = (mod 4) with D/4 = 2, 3 (mod 4) square free), let L{s, D) denote 
the Dirichlet L-series defined by 

m=l ^ ' P N i / 
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3.1. Real case. Suppose F = M, thus O = Z, do = 1 and r = 1. In this subsection we shall 
use the terminology 'quadratic form' instead M-hermitian form. Every p- local density Sp{Q, k) 
has been computed by Yang |Ya98] for every integral nondegenerate quadratic form Q and any 
/c G Z. We will use his results several times, as in the following lemma. 



Lemma 3.1. Let Q be a quadratic form as in (3.3), let k G Z, and let p be a positive prime 

m — 2 

number such that p\2 det((5) . Write k = p'^i with p\ i. Set q = p 2^ and 

'(-!)¥ dct(Q)\ 



P J 

m — 1 

(-1)^^ dct(Q)^ 



if m is even, 
if m is odd. 



Then 



(3.5) 5piQ,k)={ 



(eg) 



C+l 



(1 — ep ^ ) 
(1 -^-{"^-1)) 



eg — 1 

qc+i _ I 

g2 - 1 



if m is even, 

if m is odd and c is odd, 



q'^ — I 

(1 — — + (1 + ep^^^)q'^ if m is odd and c is even. 



g2- 1 



(-1)^ det(Q) 
P 

m — 1 



p- 



if m is even, 



1 + (-l)^dct(Q)^ ^^^^ 



In particular, if p\k, then 
(3.6) (5p(Q, k) 



Proof. By assumption, p ^ 2 and p \ ai for every z, thus /i = and = for all i, in the 
notation of |Ya98j . hence L(j, 1) = [l,m] for j even and L(j, 1) = for j odd. If m is even, 
|Ya98l Theorem 3.1] implies 

c 

(5p(Q, A;) = 1 + (1 - p-i) ^'1' - V 

c+l 

= l_(eg)^+i + (l_p-i)^(eg)^- 



(eg) 



c+l 



eg 



1-ep 2) 



The proof for even n is similar. 

The cases when p|2det(Q) can also be computed from |Ya98] . 



□ 



Theorem 3.2. Let Q be a quadratic form as in (3.3) and let /c G Z. If m is odd, write 
A;det((5) = c% (c, b E Z) with b square free and set 



D 



{-l)^b zf{-l)'^b = l (mod 4), 
4(— 1)^^6 otherwise. 



Then 
(3.7) 



C(m-l) 



n 



6p{Q, k) 



P\D 



„-(»n-l) 

p|2fcdct(Q) 



p- 
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If m is even, write det{Q) = c b (c,b eZ) with b square free and set 

D 



[-l)^b zf{-l)fb=l (mod 4) 
4(— 1) 2 b otherwise. 

Then 

Sp{Q, k) 



(3-8) W.*) = 77;-7^ n 



L(m/2,D) J-l 

^ I ' ' p|2fcdet(Q) ~ 



P' 



The proof is very similar to the proof of |RT97t Thm. 12]. It follows by decomposing k) 
as np|2fcdet(Q) ^p(<5, k) n;^2fcdet(Q) ^) ^nd by applying Lemma[ 
in the complex case. 



3.1 



We will show the details 



3.2. Complex case. Suppose now F = C, thus r = 2 and O = 7j[uj] where u = a/cIo/2 if 
do = (mod 4) or = (1 + \/d^)/2 if do = 1 (mod 4). 

In this subsection we shall use the terminology 'hermitian form' instead of C-hermitian form, 
and we denote the form by H (instead Q). Actually, in the next lemma we introduce a quadratic 
form Q associated to H. 



Lemma 3.3. Let H be a hermitian form over O as in (3.3) of signature (mi, 7722). Set 
(3.9) Qh - 



diag (ai, -oi^, . . . ,a^, -a^^) if do = (mod 4), 



diag ( ( ^J'JLu ... ;-/2 aJC'L,. ) ) ^fdo^l (mod 4). 



ai/2 ai(l-do)/4 j ' ' ' yam/2 a,nil-do)/4 ^ 

Then Qh is a quadratic form of signature (2mi,2m2) (of rank 2m) and 

S^{H,k) = 5fiQH,k) 

for every k E Z, where the upper index ¥ indicates the fields of definition of the forms. 
Proof. We have 

H[z] = ai\zi\^ ^ ha„+i|^„+iP, for z e 

Writing zj = X2j-i + X2j uj with X2j-i,X2j € Z for each zj G O, it follows that H[z] induces the 



quadratic forms in (3.9) for x G Hence 

_ i^{ze{0/p^Or:H[z]^k modf} 

# {x G iZ/p^zf"" : Qh[x] = k mod p^ 

since a representative set of O /pPO is {bi + b2UJ : bi E Z, < bi, 62 < p-'}. □ 

Remark 3.4. When p ^ 2, the matrices ^-^^2 {i-do)/'i) -^0/4) are equivalent over the 

p-adic integer numbers Zp since 

1 -l/2Vf 1 1/2 \ (1 -l/2\ ^ A \ 

1 J 1 1/2 {l-do)/4j\0 1 J I -rfo/4i- 



Consequently, both quadratic forms given in (3.9) are equivalent over Zp if p 7^ 2, thus their 



p-local densities representing any fixed k coincide. 
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Lemma 3.5. Let H be a hermitian form over O as in (|3.3|), let k E Z and let p be a positive 

number sv 
q = Then 



prime number such that p f 2da det{H) . Write k = p'^i with p ] £ and set e = ( ^ ) CLi^d 



(3.10) 



6,iH, -k) = (1 - ep-™) 



c+l 



eg — 1 



In particular, if p\ k, then 



p 



Proof. From Lemma 3.3 we see that 6p{H, k) = 6^(Qh, k) where is given by (3.9). Moreover, 
since p 7^ 2, Remark 3.4 tell us that Sp{H, k) = S^{Q'ff, k) where 



Q'h = diag (ai, -oi 



do 



4 ; • • • ; (^mi O'm 4 



)■ 



We now compute 5^ (Qj^, k) by applying |Ya98] . In its notation, since p \ do and p \ for 
all we have /j = for all 1 < i < 2m, 62^-1 = aj and = —ajdo/A for I < j < m, hence 
L{j, 1) = [1, 2m] for j odd and L{j, 1) = for j even. In particular, 1) is even for any j. 

eK Theorem 3.1 in |Ya98j now yields 



Furthermore, d{j) = —j{m — 1) and v{j) = y-y 
that 



<5(g, fc) = 1 + (1 - p-') Y,{tqy - {eqr^'p-\ 



The rest of the proof is straightforward. 



□ 



Theorem 3.6. Let H be a hermitian form as in (3.3) overO and let A; G Z. Write k = ± YlpP^'' 
with Tp e M U {0} and Tp = for almost all p. 
If m is odd, then 



(3.12) 5{H,k) 



Lim, d, 



n 



Sp(H, k) 



n 



Pi. 1 - ( T ) P-^""-'^ 

p\2dQ dct(H) 



P 



If m is even, then 



(3.13) 



6{H, k) 



n SpjH, k) yr 



1 — P 



-(Tp+l)(m-l) 



C(m) i — p 

' p\2dodet(H) ^ v\k 

p\2dQ dot(fr) 



P 



-(m-1) 



Proof. From (3.11) we have 



5{H, k) 



n 5p{H,k) n 



p\2kdo det(H) 



p\2kdo dot(H) 
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Recall that (^^^ = if P t '^o- Thus, if m is even, we obtain that 



5,{H,k)= n n (i-p^™) 

p\2kdodet{H) p\2kdo det{H) 

1 ^ 5p{H,k) 



c(itl) n 1 — »- 

^ p|2fcdc)det{_H') ^ 



1 



n SpjH, k) yr a 



C(m) 



Sp{H, k) 



p\2do dct{H) 



p\2d(j dct(H) 



Hence, (3.13) follows from (3.10). 



On the other hand, if n is even, we have that 

sp{H,k)= n WA:) n (^-(t)^"™) 

p\2kdodet{H) p\2kdodet(H) ^ V / / 

1 T-r 5p(i/,A;) -pr 



n 



L(n + l,do) 1 Mo 



n 



p\2dQ dot(H) 



p- 



Finally, (3.12) follows from (3.10). 



□ 



(4.1) 



4. Numerical approximation of the error 

The goal of this section is to approximate the behavior of 

Nt{Q,-k) 



<^{t) = ^t{Q,-k) 



f2p 



CiQ,-k) 



for t < T, where T = 10^, lO'' if F = M, C respectively. Here and subsequently, we consider 
G N and Q = an F-hermitian form as in (1.8) where a G N and A is diagonal with 

positive integer entries ai, . . . , a„, that is 



(4.2) 



Q 



\ 



-a) 



This section is divided into three parts: calculation of C{Q,—k), calculation of Nt{Q,—k) 
for every t < T and finally, approximation of the behavior of ^E'(t). 



4.1. Calculation of C(Q, —fc). By writing the main coefficient (1.11 ) as C((5, —A;) = C'{Q)6{Q, 
C'{Q) depend only on the F-hermitian form Q. If Q is as in (4.2) then 

2 7r"/2 



(4.3) 



n-2 

a 2 



C'iQ) = < 



(^-l)r(f) (ai...a„)i 

2^+1 ^n+l ^n-l 

L n! Idol 2 ai . . .an 



if F 



if F = C, 



since r = dimR(F), 2p = n - 1, 2n if F = M, C respectively, and vol(5™~^) = 27r'"/Vr(m/2). 

Furthermore, the local density 6{Q, —k) has been calculated in Section [3| up to finitely 
many p-local densities 6p{Q,—k). However, all these terms are computed in |Ya98] . In the 
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complex case we have to use Lemma 3.3 We usually use the advanced quadratic forms library 
in Sage [Sage], programmed by J. Hanke and A. Haensch. 

27rV3 



For example, for Q = Ia^i and k = 1 considered in Example 



by (4.3) and, by Theorem 3.2 



1.4 



L(2,4) <52(/4,i,-l) vrV8 5/1 



we have that C'{l4^i] 
15 



Hence C(J4, 



C(4) 

1) = 5. In Example 



1 



7rV90 15/24 2 



■7r 



-2 



1.5| we take H 
that C'iQ) = (27r)V(5! 3^) by Q and 6{h,i 



h 1 over Z 



and k = 1, then we have 



_ 1) = 1/C(6) . 52(/5,l, -l)/(l-2-^) . 53(/5,l, 

3^6) by ( [3l2| . We check at once that C(J5,i, -1) = 18. 

4.2. Calculation of Nt{Q, —k). We consider the functions FA,oiTn) '■= ^ C^" : = 1} 
and Go{fn) := G (9 : |zp = m}. It follows by definition ( 1.9[ ) that 



(4.4) 



7Vi(Q, -fc) = J] FA{am - k) Go{m) 



m=l 



In the real case, (4.4) reduces as Nt{Q, —k) = 2 X]m=i FA^arn? — k). 

There exist formulas for FA,o{fn) cind Go{fn) for some choices of A and O, that allow us 
to have fast algorithms to compute Nt{Q, —k) for large t. For example, we have used Jacobi's 
four-square formula 

d\m,A]d 



in Example 1.4 The other formulas for ^^(m) that we use here can find in |AAW08] . |AAW10] 
and jOtfl] . 

4.3. Estimation of the decay order of "^tiQ, —k). We have computed in the last two sub- 
sections the necessary elements to obtain the function given in (4.1). This function goes 
to zero when t — )■ +oo. Our purpose is to estimate the order of decay of this function. The 
method used here is new as far as we know. 

We consider the set of points in the plane given by 

S = {(m,^(m)) G |1,T] X M : ^(m) > ^(/) V/ > m} . 

Now, by the "least square fitting-power law" method, we approximate these points by a curve 
y{t) = Bt" with i? > and a < 0. This method works as follows: if {(xj, ?/i)}i<j</i are points 
with positive entries, then B and a are given by 

h h h 

h^{\ogXi\ogyi) - J^(logXi) ^(log?/i) 



i=l 



i=l 



i=l 



B 



h^[\ogXi)'^ - (^^log 

i=l i=l 
\i=l i=l / 



Xi 



Finally, our estimate of the decay order of \l/(t) is 0(t'^). 

For instance, in Example 1.4 we obtain the curve y{t) = 13t~°'^^ and in Example 1.5 the 
curve y{t) = 56 1"^ '^^. 
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5. Results and conclusions 



We conclude the paper by showing the experimental results in several tables, together with 



some conclusions. Recall that formula (1.10) proved in |Lal2j implies that "^{t) = "^tiQ, —k) ^ 
f^-^p ^Jigj-g J- ig defined by (1.3). Hence limt_>.oo "^(t) = since r < 2p. 



In Corollary |2.4| and Remark |2.7| we obtain that 

-1/3 ifF 
r-2p< T := <( -1/2 ifF 
-1 ifF 



M and n = 2, 
M and n > 3, 
C and n > 2. 



Recall that, in the complex case, we are assuming the extended result of Theorem 2.6 Similarly, 



from (1.7), by assuming that condition (1.6) holds (e.g. there are no exceptional eigenvalues) 
then 

2 

-1 + , 

T — 2p = \l :- 

-2 + 



ifF 



n + 1 



+ e 



if F = C, 



for each e > 0. On the other hand, we have estimated the decay order of "^tiQ, —k) by a curve 
y{t) = Bf^ (see Subsection |4.3 ). Now, our purpose is to compare the number a with T and VL. 
In Tables [Tl [21 13] and |4] we consider the M-hermitian forms (or quadratic forms) 



Q 



with 1 < a < 15, 



for n = 2, 4, 6, 8 respectively (/„ denotes the n x n-identity matrix). Similarly, Tables [sj [6| [7| 
8] consider the C-hermitian forms (or just hermitian forms) H = diag(J„, —a) with 1 < a < 15, 
over the ring O = Z[-\/— 3] for n = 2,3,4,5 respectively. In each table is fixed k = 1, and 
T = lOMf F = M and T = 10^ if F = C. 

We abbreviate the results from Tables [THE] as follows: 

F n T„ J- Qn,¥ cr in Table 



(5.1) 



2 -0.33 -0.33 

4 -0.50 -0.60 

6 -0.50 -0.71 

8 -0.50 -0.78 



[-0.68 
[-0.98 
-0.999 
-1.002 



C 



2 
3 
4 
5 



-1 
-1 
-1 
-1 



-1.33 
-1.50 
-1.60 
-1.67 



[-1.40 
[-1.37 
[-1.37 
[-1.37 



-0.50] 
-0.93] 
-0.991] 
-0.9976] 



-1.33] 
-1.32] 
-1.33] 
-1.34] 



We divide the conclusions according to a is larger or not than Q, or equivalently, the numerical 
data gives evidences of exceptional spectrum or of a new lower bound for Ai. 

5.1. Evidences of exceptional spectrum. We see from Tables [6]-[8] that a > Qn,c for any 

n = 3,4,5, hence these approximations say that (1.6) should not hold, in particular, there are 



1.2 



exceptional eigenvalues in rQ\H^. This numerical data supports Conjecture 

Let V = Vn,c = 2n — 1 and u = Un,c = '^n? /{n + 1)^. Hence, Ai > is the extended lower 
bound proved in Theorem 2.6 and Ai > w is the condition (1.6). When n = 2, since a < f22,c 



then Ai > u. For 3 < n < 5, a lies (approximately) in [—1.4,-1.33] by (5.1). Hence, these 
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approximations tell us that r — 2p = a/p^ — Ai — p lies in the same interval, or equivalently, Ai 
lies in the intervals showed as follows: 



,2 



n V oj Ai in p 

~2 3 3^55 [3.55, 4] T 

(5.2) 3 5 6.75 [6.22, 6.44] 9 

4 7 10.24 [8.88, 9.24] 16 

5 9 13.88 [11.55, 12.04] 25 








V 


— ♦ ♦ 

UJ \i p 





V 


e— • 

XiUJ 


P 





• e— 

V Ai 

• e — • — 


UJ 


P 



V Xi UJ 



5.2. Evidences on a lower bound for Ai. We now restrict our attention when a is smaller 
than Vt and T, namely, the cases F = M with n = 2,4,6,8 and F = C with n = 2. As we 
explained in Introduction, these data tell us that the error term of ( 1.10 ) is better than the error 
term expected by the lattice point theorem, even by assuming the nonexistence of exceptional 
eigenvalues in FgyHg. Therefore, these computations are numerical evidence that supports 
Conjecture L3 in these cases. 

Some "noise" is visible in the computation of = 2. Actually, the amplitude of the values 
taken by a is dramatically smaller as n grows. This noise makes us doubt about how good is 
the estimation of the error term in (1.10) when n = 2. However, we do not conjecture 

anything in this case [n = 2) since the lower bound obtained in [KS03] gives the best error 
term in (1.10) that we can obtain form |Lal2j . In other words, when F = M and n = 2, 
we are showing evidences for the lower bound (1.6), namely Ai > 2/9 = 0.22---, which was 
already proved by Kim and Sarnak |KS03] . namely Ai > 975/4096 = 0.238 ■ ■ ■ . Perhaps, the 
complications that emerge in this case (n = 2) are connected with the same complexities of the 
Gauss circle problem. 

We actually consider some other F-hermitian forms that are not included here. They are 
diag(l, 1, 2, 2, —a), diag(l, 1, 1, 4, —a), diag(l, 1, 3, 3, —a), diag(l, 2, 2, 4, —a), diag(l, 4, 4, 4, —a), 
diag(l, 1, 4, 4, — a), diag(l, 1, 1, 1, 2, 2, — a) and diag(l, 1, 2, 2, 2, 2, — a) when F = M and also 
diag(/„, —a) over for n = 2, 3, 4 when F = C, with 1 < a < 15. We have used formulas 

for FA{m) from |AALW07] and |AAW08j . The numerical results in these cases are very similar 
those exposed here. 
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Q = diag(/2, -a) 
F = M, n = 2, T 



■2,: 



-0.33, 

a C(Q,-1) 



1 < a < 15, 

= 10^ 

= —0.33. 



1 


2 


2 


AV2 


3 




4 





5 


|V5 


6 


4x/6 


7 





8 





9 


4 


10 




11 




12 





13 




14 




15 






0.000492 
0.001647 
0.000691 

0.000427 
0.002946 



0.001732 
0.001229 
0.002404 

0.001848 
0.003855 



-0.51 
-0.56 
-0.65 

-0.68 
-0.60 



-0.60 
-0.57 
-0.57 

-0.50 
-0.50 



Q = diag(/4, —a) 
F = R, n = 4, T 



T 



= -0.5, 

a c(g,-i) 



1 < a < 15, 

= lo^ 



-0.6. 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



5 

8^2 
12^3 

32 
f^/5 

fVe 

32\/2 
36 

280 /Jo 
13 ^ 

f v/3 

364 /T^ 
17 V J-<3 
112 . 

5 
360 

13 



14 
l5 



0.000004 
0.000015 
0.000032 
0.000039 
0.000030 
0.000059 
0.000082 
0.000060 
0.000008 
0.000072 
0.000164 
0.000101 
0.000073 
0.000136 
0.000200 



-0.93 
-0.97 
-0.98 
-0.95 
-0.98 
-0.98 
-0.98 
-0.97 
-0.97 
-0.98 
-0.97 
-0.97 
-0.97 
-0.97 
-0.98 



Table 3. 



Table 4. 



Q = diag(/6, -a) 

F = M, n = 6, T 

T6,M = -0.50, 



1 < a < 15, 

= lo^ 

^6,R = -0.71. 



a C{Q,-1) *(T) 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 



"27" 
4 

18V2 

945 

26 V <5 
120 

360^2 

15309 

26 

3950 /Tn 
21 V -LU 

9801 /3J 
38 ^ 

689013 
2198 

6468 /;L4 

19 ^ 
5550 . /iT 



13 



0.000016 

0.000048 
0.000135 
0.000327 
0.000400 
0.000450 
0.000832 
0.001283 
0.001346 
0.001198 
0.001887 
0.002874 
0.002703 
0.002454 
0.003827 



-0.992 

-0.998 
-0.998 
-0.991 
-0.998 
-0.999 
-0.998 
-0.999 
-0.996 
-0.998 
-0.997 
-0.998 
-0.997 
-0.998 
-0.998 



Q = diag(/8, -a) 

F = M, n = 8, T 

^8,M — —0.5, 



a 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



1 < a < 15, 

= lo^ 

f^s.K = —0.78. 



C(Q,-1) *(T) 



"585" 
136 

46575 . 
^120'^'^ 

#'x/5 

281880 /c 
697 
11627700 /j 

20417 
™^ 

2496825 



697 
7885000 
5321 
214923225 
124457 
3179520 AT 
697 V O 
43006275 

14281 
68682320 

20417 
815568750 

218161 



13 

Ti 

l5 



0.000018 

0.000127 
0.000377 
0.001031 
0.002218 
0.003478 
0.004801 
0.008191 
0.013384 
0.015999 
0.018761 
0.027650 
0.039000 
0.044185 
0.047254 



-0.9976 

-0.9999 
-1.0000 
-0.9991 
-1.0001 
-1.0001 
-1.0001 
-1.0002 
-0.9996 
-1.0001 
-0.9999 
-1.0001 
-1.0001 
-1.0002 
-1.0001 
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Table 6. 



Q = diag(/2, —a) 1 < a < 15, 
F = C, = n = 2, T 

'^2,c = —1, ^2,c = —1.33. 



10^ 



Q = diag(/3, —a) 1 < a < 15, 
F = C, = Z[v^], n = 3, T 

Ta.c = -1, ^3,c = -1-5. 



10^ 



Q 


C{Q-V) 


^\J- } 


(J 


(X 


1 


18 


000039 


-1.33 


1 


2 


48 


OOOOfil 


— 1 34 





3 


108 


000102 


— 1 36 


3 


4 


48 


000147 


-1.37 


4 


5 


1200 


0.000050 


-1.35 


5 


6 


144 


0.000255 


-1.34 


6 


7 


2352 
19 


0.000196 


-1.40 


7 


8 


192 


0.000243 


-1.35 


8 


9 


324 


0.000100 


-1.37 


9 


10 


800 

14^20 
37 


0.000296 


-1.35 


10 


11 


0.000686 


-1.36 


11 


12 


288 


0.000345 


-1.34 


12 


13 


14196 


0.000571 


-1.36 


13 


14 


0.000388 


-1.37 


14 


15 


3J80 

7 


0.000407 


-1.37 


15 



30 

576 
5 

324 

576 

7875 
7f?6 
58^53 

2916 

47250 
88^46 

3??e4 



1206153 

7(fef36 

21^^25 
26 





(J 


000085 


-1.35 


000349 


-1.32 


000884 


-1.35 


001443 


-1.35 


n 0091 QQ 


— 1 

-L . OU 


0.003006 


-1.36 


0.004731 


-1.36 


0.005813 


-1.34 


0.007188 


-1.36 


0.008708 


-1.35 


0.010820 


-1.36 


0.011752 


-1.37 


0.014859 


-1.35 


0.017702 


-1.34 


0.022768 


-1.36 



Table 7. 



Table 8. 



Q = diag(/4, —a) 1 < a < 15, 
F = C, O^Z[y^], n^4:, T 

Tac — —1, ^4,C = —1.6. 



10^ 



Q = diag(/5, —a) 1 < a < 15, 

F = C, O = n = 5, T = 10^ 

T5,c = -1, ^^5,c = -1.67. 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



C{Q,-1] 



30 

1920 
11 

648 

19200 



156^)1 



00 



5fiio 

2881^000 



12 



go 



17496 

156000000 

5731 
428688480 



14 



2039^^5400 



18. 



QOO 



^mo 



421 



521 





a 


0.000052 


-1. 


.36 


0.000468 


-1 


.33 


0.001536 


-1, 


.35 


0.003596 


-1, 


.35 


0.006766 


-1, 


.35 


0.012691 


-1 


.36 


0.018688 


-1 


.36 


0.029712 


-1, 


.36 


0.039107 


-1, 


.36 


0.057215 


-1, 


.36 


0.076104 


-1 


.36 


0.100624 


-1, 


.36 


0.111032 


-1, 


.37 


0.158689 


-1, 


.37 


0.189703 


-1 


.36 



a C{Q,-1) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



18 

4224 
17^)6 



488I375 



2otOg64 

141^2^221 



1« 



44 



I4M76 

2865'!?0000 



1291 



'^2826 



32faf24 
34464^0139 



82i 



632 



237^8^(^50 
2821 





a 


0.000059 


-1. 


.36 


0.000920 


-1 


.34 


0.004840 


-1, 


.35 


0.014810 


-1, 


.35 


0.036694 


-1, 


.36 


0.075213 


-1 


.36 


0.142486 


-1 


.37 


0.234356 


-1, 


.36 


0.387823 


-1, 


.36 


0.577332 


-1, 


.37 


0.866876 


-1 


.37 


1.201558 


-1, 


.36 


1.680102 


-1, 


.37 


2.199171 


-1, 


.37 


3.014865 


-1 


.36 
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